N84 2V310 




MHTHOD TO ESTIMATE DRAG COEFFICIENT AT THE 
AIR/ICE INTERFACE OVER DRIFTING OPEN PACK ICE 
FROM REMOTELY SENSED DATA 


Uri Feldman 

Department of Geography, Bar-Ilan University 
Ramat-Gan 52100, Israel 


ABSTRACT 


A knowledge in ne'’r real time, of the surface drag coefficient for drifting pack ice 
is vital for predicting its motions. And since this is not routeinly available from 
measurements it must be replaced by estimates. Hence, a method for estimating this 
variable, as well as the drag coefficient at the water/ice interface and the ice 
thickness, for drifting open pack ice was developed. These estimates were derived 
from three-day sequences of Landsat-1 MSS images and surface weather charts and 
from the observed minima and maxima of these variables. "Hie method was tested with 
four data sets in the southeastern Beaufort sea. Acceptable results were obtained 
for three data sets. Routine application of the method depends on the availability 
of data from an all-weather air or spaceborne remote sensing system, producing images 
with high geometric fidelity and high resolution. 


1. INTRODUCTION 

The study of drifting pack ice depends on the knowledge of numerous parameters 
related to its motion (feldman and Howarth, 1979). And since data vital for deter- 
mining these parameters are not routinely available for the polar oceans they must 
be estimated (Feldman et al., 1979, 1982). The purpose of this study is to present 
a method for estimating ice thickness and drag coefficients at the air/ice and the 
water/ice interfaces for groups of detached ice floes. To this end, a reduced form 
of the general equation of motion for drifting pack ice was employed, assuming wind 
stress, water drag and Coriolis force to be at equilibrium. 

Data were obtained from three sources: 

(1) Pack ice speed and direction of motion, which were measured fr: 'hree-day 

sequences of sidelapping Landsat-1 mullispectral scanner (MSS) :,es. 

(2) Surface wind speed and direction, at 10m above the ice and air density, which 
were obtained from three-day sequences of surface weather charts, following 
the work by Feldman et al. (1979, 1981). 

(3) Minimum and maximum values of pack ice thickness and drag coefficients of 
its surface and subsurface, which were obtained from observations previously 
conducted in the Arctic ocean. 

The method developed in this study is based upon employing three latios between 
pairs of the unknown parameters, which could be calculated from the pack ice velo- 
city vector and the wind field data, in conjunction with the known minima and maxima 
of the same parameters. The procedure to obtain the estimates was tested by four 
groups of open pack ice, consisting of three or four data sets within each group, 
which drifted in the Beaufort sea during October 1973, July and August 1974. 
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Acceptability of the results was determined from the corresponding values of the 
cross isobar angle. 


2 . PHYS ICAL BACKGROUND 

The general equation of motion for a unit area of drifting pack ice is given 
by Campbell (1965} as 

T +T +C+P+T=p. h dV/dt (1) 

a w ice 

where t is the horizontal air stress at the air/ice interface, t is the hori- 
a ^ w 

zontal water stress at the water/ice interface, C is the horizontal Coriolis 

— V — V- 

deflecting force, P is the horizontal marine pressure gradient force, I is the 

horizontal internal ice stress, p. is the ice density, h is ice thickness and 
^ ice 

dV/dt is the horizontal ice acceleration. Thorndike (1973) has shown that the 
horizontal acceleration term is usually much smaller than all other terms in 
equation (1). Therefore, a steady state drift can be assumed (Nansen, 1902; 
McPhee, 1982) and equation (1) can be rewritten as 

t+t+C + P+ T = 0 (2) 

aw 


If this equation is applied to drifting open pack ice, consisting of detached 
ice floes, where internal ice stress cannot be transmitted among the floes 
(Hibler, 1979; McPhee, 1980) then equation (2) can be reduced to 

T + T + C + P = C (3) 

aw 


Using data reported by Newton (1973, p.23) it can be shown that mean speed 
of ocean currents in the Beaufort sea is less than 2m day^. In addition, Hibler 
and Tucker (1979) stated that geostrophic currents and ocean tilt have a negligible 
effect on short term, weekly drifts. Hence, it may be concluded that in the area 
of study motions of drifting open pack ice can be determined from a simple steady 
state equation, written as 


T„ + T, + C = 0 
a w 


( 4 ) 


Nansen (1902), Sverdrup (1928), Shuleiken (1938), Fel’zenbaum (1958), 

Campbell (1965), Thorndike (1973), Neralla et al. (1980), Feldman et al. (1981) 
and McPhee (1982) applied equation (4) in their studies. 


Under conditirns of neutral equilibrium within the atmospheric boundary layer 


urn 


a 




( 5 ) 


where p is the air der ity, C, the drag coefficient at the air/ice interface 
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and U the horizontal surface wind speed at 10 m above the ice surface. Under 
similar conditions beneath the ice (Johannessen, 1970) 

T = p v2 (6) 

w w d 

w 

where is the ocean v/ater density, is the drag coefficient at the water/ice 
interface and V is the speed of the centre of gravity of a drifting group of ice 
floes. The Coriolis force may be derived from 

C = p. f h V (7) 

ice 

where h is the ice thickness and f(=2cj sin(J)) is the Coriolis parameter, w 
(=7.292 10"^ s'^) is the Earth’s angular speed and (}> is latitude. 

-V -y 

Resolving the x and y components of t , and C from equations (5), (6) and 
(7), allows equation (4) to be written as 

p d cos Ay = p Cj V2 (8) 

’^a d ’ ^w d 


^d "" '^ice ^ ^ 

where Ayi the angle of sea ice deflectio.n (Feldman et al . , 1981) is defined as 

Ay = e. - 9 (10) 

' ice u ^ ^ 

and 0. and 9 are the directions of motion of the pack ice and the surface wind 
ice u 

respectively. 

The cross isobar angle, A9 is defined as 

A9 = 0^ - 0^ (11) 

where 9 q is the geostrophic wind direction. A0 may be obtained from the dif- 
ference between equations (11) and (10), written as 

“ = «G - ®ice 


3. PACK ICE VELOCITY FROM LANDSAT MSS IMAGES 

A number of techniques have been used to calculate the velocity of drifting 
pack ice from sequential Landsat MSS imagery (Crowder et al , , 1973; Hibler et al., 
1974; Wendler and Jayaweera, 1974; Nye and Thomas, 1974; Nye, '•975 and Sobczak, 
1977). The orbits of Landsat converge in high latitude thereby producing sequences 
of four sidelapping images over the Beaufort sea. In this study, velocities of 
four groups of drifting ice floes were calculated over three-day sequences. A 
group could consist of any number of single ice floes in close proximity, but in 


ji-. 



these cases they ranged from 3 to 26. To determine the velocity it is necessary 
to measure the co-ordinates of each ice floe in a group and to know the exact 
mean time of imaging, tj, on each day. The co-ordinates, x^; yj, of the ice 
floes were measured with a digitizer and were related to an origin and several 
control points located on land. The exact mean scanning time, ti, was determined 
from the orbital information. 

The area, Ai, of each floe was calculated from 1:250,000 scale enlargements 
of the images, using the dot grid method. The co-ordinates of the estimated 
centre of gravity of each group, X ; Y , were calculated from 

V " ^ V " ^ 

This procedure eliminated effects due to collisions which could occur within a 
group while in motion. 

The component mean velocities of drifting centres of gravity, V^; Vy, were 
calculated for the intermediate scanning time, ^"^12* ^27’ ^34 ^ from 


\ CIS) and Vy = AYg^/At ( 16 ) 

where At is the time increment between sequential paths of Landsat-1 and AX or 
AYg^ are the component distance increments during At. 


4. RATIOS AND MEAN RATIOS BETWEEN ICE PARAMETERS 

a w 

The first step towards estimating h, Ch d Cj consisted of determining the 
ratios M, N and B and their means M, N an'’ • lues of M and N, defined as 

M = h/C^ (17) and N = C^/C^ (18) 

were calculated from available data, for the four groups, at by rewriting 

equations (17) and (18) from equations (9) and (8) respectively, as 


M = (p^ U2 sinAY)/(Picg ^ 

(19) 

N = (p^ U2 cosAy)/(p^ V2) 

(20) 


The variable U was obtained from the geostrophic wind speed G by the formula 

U » 0.54G + 1.68 (21) 

adopted by Feldman et al. (1979) from Basse's (1974a, 1974b) work. G and were 
derived from surface weather charts. 
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Ay was replaced by values of the mean angle of ice deflection Av, defined 
and calculated by Feldman et al. (1981) in the area of study for V was 

calculated from equations (15) and (16). Finally, the constants <p (required for 

calculating f in f = 2 w sin (P) and were replaced by 4i - 70“, = 1.08 

10^ kg m"^ and p. = 0.91 10^ kg m“^ respectively. 

X C 6 

The means M and N were calculated for the four data sets at t ^ 5 t^,, at t.,- 
5 t^^, at t ^2 5 ^34 tj 2 » ^2.8 ^ ^34' of B, defined as 

B = h/Cj (22) were calculated from: B = ,M/.\' (23) 

by replacing M and N in equation (23) with M and N. Values of B, M and \ are 
presented in Table 1. 


5. ACCEPTABILITY OF RESULTS FOR B, M AND N 

_ _ Resul_ts presented in Table 1 indicate that two fully acceptable ratio sets of 
B, M and N are "vailable for cycle 26a, as well as one for cycle 26b, one for cycle 
41 and none for cycle 43. And although the acceptable ratio sets obtained are 
sufficient for deriving the estimates of h, Cj and C^, for three out of the four 
cycles tested, it is evidently necessary to account for the frequent occurrence 
of ratio sets or ratios which were either not available or rejected in Table 1. 

^ Availability of Landsat-1 MSS images 

Ratio sets numbers 4.2, 4.3 and 4.4 could not be determijied because the 
Landsat-1 image, required for obtaining the variables V and for cycle 43 at 

t^, was not available. This might occur in cases where a dense cloud cover pre- 
vents identification of ice floes on an image or in cases where floes drift outside 
the area covered by the corresponding image. 

5.2 Ratios rejected by A? 

Ratio sets number 2.1, 2.2, 2.3, 3.2, 3.3 and 3.4 were rejected where values 

of A0, calculated from equation (12) were either less than 0* or greater than 60®. 
This range of A0 was chosen as the criterion of acceptability for the ratio sets, 
firstly because corresponding observed data, required for calculating M and N from 
equations (19) and (20) respectively, were not available for comparison and 
secondly, because this range, which was determined from observations, is relatively 
small. The limits of Jd (0‘-and 60®) were determined from studies by Gordon (1952), 
Reynolds (1956), Aagaard (1969), Masse (1974a, 1974b) pnd Lavrov (1974) conducted 
over sea surfaces, which have drag coefficients similar to those over pack ice 
(Roll, 1965) and from studies by Smith et al. (1970), Banke and Smith (1973), 

Feldman et al. (1979) and Albright (1980) conducted over pack ice. 

Data of V and used in this study, are considered to be highly accurate, 

because Landsat-1 MSS images are nearly free of distortions and because the tech- 
nique employed to obtain these data produces accurate results. Hence, rejection 
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Table 1; Means of the Ratios, B, M and N 
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(— ) Ratios not )?iven: Where data on pack ice velocity are not available. 



of a ratio may result either from errors in the variables Ay, fi,; and fequations 
19 and 20) and/or from errors in ^ and 0q (equation 12). And since variations 
in Ay (Feldman et al., 1981) and Oa ^^^e relatively small, .and U is linearly re- 
lated to G (equation 21)_^ it follows that rejection of a ratio set is mainly due 
to errors in and/or 0 q. Employment of interpolated values of G and which 
were needed to replace gaps in the data sequences obtained from the surface weather 
charts, could be the main source of error. 


6. MINIMA AND MAXIMA OF RATIOS BETWEEN OBSERVED ICE PAR/\METERS 

3. W 

The second step towards estimating h, C^j and Cj consisted of determining the 
observed minima i.nd maxima Bq, Mq and N^ from the observed minima and maxima ho, 

Cjo atid C^o* which had previously been measured in the Artie ocean by other investi- 
gators. The ranges of hp, and were summarized by Feldman et al. (1981) as 


0.00 = 


h min 
0 



h 

0 

4 


h 

0 max 

= 3 . 00 m 


(24) 

0.95 = 

103 

C'j • 

do min 

4 

103 

C^ 

do 

4 

103 

do max 

= 4.00 

and 

(25) 

3.32 = 

103 

c'" 

do min 

4 

103 

c” 

do 

4 

103 

do max 

= 57.17 


(26) 


Minima and maxima of Bq, M and N^ we^e determined by replacing h, Cj and 
Cj in equations (22), (17) and (18) with the observed minima and maxima or these 
parameters, given in equations (24), (25) and (26). These were 

0.00 4 h i 903.61 m (27) 

o 

0.00 ^ M 4 : 3157.89 m and (28) 

o 

0.83 ^ N « 60.18 (29) 

o 


7. ESTIMATING ICE THICKNESS AND DRAG COEFFICIENTS 

The final stage of estimating ice thickness and drag coefficients in the area 
during the period of study consisted of rewriting equations (22) and (17) for h, 
using the observed minima Cjq , Cjq andjh^ from equations (26), (25) and (24), 
respectively, and the calculated means B and M, written as 


h = Cjo B > 3.32 10-3 B „ (30) 

h » Cjo M > 0.95 10-3 M m and (31) 

h « h > 0.00 m (32) 

0 

h > max. (3.32 10-3 f, 0.95 10' 3 M, 0.00} » m (33) 
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hence 


4 



where the lower limit of h, was given by the maximum (’max.) among the three 

values of h in equation (53). In the same way hyi^, the upper limit of h was 
given by the minimum (min.) among the three values of h in equation (34). 

hvf min. {57.17 10'3 f, 4.00 10'3 M, .5.00} = m (34) 

The range of h was defined as 

Ranges_of h were determined for each data set from equations (33) and (34) with 
B and M from Table 1. Results are presented in Table 2. 

a w 

Ranges of Cj and Cj were derived from (17) and (22) after replacing M and B 
by their means and h by its lower and upper limit, written as 

h^i/M 4 4 h^j^/M (36) and h^^^^/B ^ c” h^j^/B (37) 

a w 

Ranges of Cj and Cj were determined for each data set from equations (36) and 
(37) respectively with B and M from Table 1 and h^j^ and h^^ from Table 2. The 
results were presented in Table 2. 


8. EVALUATING THE ESTIMATES OF h, AND cj[ 

When the calculated values B, M and N (Table 1) are either less than the 
corresponding minimum or greater than the corresponding maximum of Bq, Mq and Nq 
( equations 27. 28 and 29 respectively), then the estimated values of the lower 
limit of h, C| and are greater than the upper lim: zf these variables. In 

these cases the results are unacceptable as estimates for h, CQ and C^. Hence, 
the results for cycle 43 in Table 2 were considered unacceptable. 

_ The best results were those obtained for cycles 26b and 41, where E, M and 
N (Table 1) are either greater than the minimum or less than the maximum of Bq, 

Mq and N© (equations 27, 28 and 29 respectively). For the same reasons the results 
<^or cycle 26a are partially acceptable (i.e., the results for ti2 tt,d and for 
^12’ ^23 ^ ^34 acceptable while those for t ^2 5 ^23 ^23 ^34 not)* 

The calculated values of the ranges of h obtained for cycles 26a and 26b do 
not agree with those measured at the nearest coastal stations (Table 2). Freeze up 
along the west coast of the Mackenzie bay occurred about 8 October 1973 (inter- 
preted from Landsat-1 image 1442-20295) and at Tuktoyaktuk, N.W.T. on 9 October 
1973 (AES, 1974). Since only 0.25- 0.50 m of ice could be formed during a period 
of 2-3 weeks (Pounder, 1965), it was suggested that values of h, measured during 

cycles 26a and 26b, provide a better approximation to the actual values than the 

calculated one. The calculated value of hLi^ for cycle 26b is less than hLL 
■* cycle 26a and provides, therefore, a better estimate fer the actual value of h. 

^ The calculated values of h for cycle 41 were within the range of values measured 

^ at Sachs Harbour, N.W.T. (1.80m) on 14 June 1974 and at Cape Parry, N.W.T. 

- (1.07 m) on 5 July 1974 (AES, 1974). 
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CONCLUSIONS 




It has been demonstrated that ice thickness and drag coefficients, associated 
with drifting open pack ice and vital for predicting its motions, which are not 
routinely available for the polar oceans, can be estimated from three-day sequences 
of satellite images and wind field data. Unfortunately, the method used could not 
be tested by a larger number of data sequence.^, since Landsat-1 has an 18 day re- 
peat cycle and data recorded on cloudy days are useless. 

Images recorded by weather satellites, w.thin the visible or the thermal 
infrared regions of the electromagnetic specurum were considered for application 
to this study. They were rejected because their images are geometrically distorted 
and have a relatively low resolution. 

It is recommended that for a routine application of this method an all weather 
air or spiaceborne remote sensing system, producing high geometric fidelity and high 
resolution images, should be employed. 
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